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ABSTRACT 

We report the discovery of planetary-mass companions to two red giants by the ongoing Penn 
State - Torun Planet Search (PTPS) conducted with the 9.2-m Hobby-Eberly Telescope. The 1.1 Mq 
KO-giant, BD-l-15 2940, has a 1.1 Mj minimum mass companion orbiting the star at a 137.5-day 
period in a 0.54 AU orbit what makes it the closest - in planet around a giant and possible subject 
of engulfment as the consequence of stellar evolution. HD 233604, a 1.5 Mq K5-giant, is orbited 
by a 6.6 Mj minimum mass planet which has a period of 192 days and a semi-major axis of only 
0.75 AU making it one of the least distant planets to a giant star. The chemical composition analysis 
of HD 233604 reveals a relatively high ^Li abundance which may be a sign of its early evolutionary 
stage or recent engulfment of another planet in the system. 

We also present independent detections of planetary- mass companions to HD 209458 and HD 88133, 
and stellar activity- induced RV variations in HD 166435, as part of the discussion of the observing 
and data analysis methods used in the PTPS project. 

Subject headings: planetary systems-stars: individual (BD+15 2940, HD 233604, HD 209458, 
HD 88133, HD 166435) 



1. INTRODUCTION 

A precision of the radial velocity (RV) technique, which 
is the most prolific ground based method of planet detec- 
tion, depends on the presence of a large number of narrow 
spectral lines in the stellar spectra. This limits its effi- 
ciency to slowly rotating, cool stars of late-F to M spec- 
tral types which also sets an upper limit of Mi, < 1.5M0 
to the masses of Main Sequence (MS) target stars that 
may be hosts to planetary companions. Fortunately, as 
these stars evolve off the MS and become red giants, they 
cool and slow down their rotation which makes them suit- 
able again for precision RV measurements. Despite the 
fact that the precision of this technique is se riously af- 
fected by stellar activity (jOueloz et al.l 120011) which is 
particularl y true in the cas e of giant stars (jCarnev et al.l 
[2003; Scti awan et al.|[2004D . it remains most efficient in 
searches for planets around A-F stars with masses as high 
as 5 Mq . This offers an excellent opportunity to improve 
our understanding of the population of planets around 
evolved stars that are significantly more massive than 
the Sun. 

Searches for planets around red giants also represent 
an essential extension of the traditional surveys, because 
each detection of a planetary system around a suffi- 
ciently evolved giant provides a snapshot of the changes 
in its dynamical configuration powered by the evolution 
of the parent star. In principle, this makes it possi- 
ble to track the evolution of planetary systems back to 
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their origins (jAdamow et al.l I2012D and to extrapolate 
it beyond the r ed giant stage ([Villaver and Livid [20071 : 
iMustill and Villaver 201^ . 

Over a dozen stars with planetary-mass companions 
have already been discovered by the Penn State- Torun 
Centre for Astronomy Planet Search (PTPS), whose pri- 
mary, long-term goal is to improve our understanding 
of the evolu tion of planetary systems around aging stars 
([Niedzielski et al.ll2007l l2009al lbl: iGettel et al.ll2012 b.al. 
Our findings are part of the growing body of discover- 
ies by s urveys such as the McDonald Observatory Planet 
Searc h (jCochran and HatzesI 119931: IHatzes and CochranI 
I1993D . Okayama Plane t Search (ISato et al.l l2003D. TauT 
enburg Planet Search (" Hatzes et al.l 12005). Lick K-giant 
Survey (Frink et al. 2001 , ESO FEROS planet search 
(jSetiawan et al. 20031). Retired A Stars and Their Com- 



panion s (iJohnson et al.l 



2007), CORALIE & HARPS 



search ([Lovis and Mavoii[2007D , Boyunsen Planet Search 
(jOmiva et al.ll201ll) . and several others. 

The existing population of planets around evolved 
intermediate-mass stars, uncovered by these surveys, 
begins to show its distinction from that of the solar 
type MS dwarfs. For example, the planet frequency- 
stellar metallicity correlat ion observed for MS dwarfs 
(|Fischer and Valeiitl I2005D appea rs to be weaker in 
the case of planets ar ound g iants (IPasauini et al.l 120071 
iHekker and AIelende3 [20071: iGhezzi et al.l I2010D . and, 
compared to dwarfs, the frequency of i ntermediate-mass 
stars with planets seems to be hi gher ([Lovis and Mavoij 
I2007t IJohnson et alll2010al 1201 ID . Moreover, the inter- 
mediate mass stars prefer to have more massive planets 
and are deficient in low mass (< 1.5 M,j) companions 
(ISato et al.ll2007l: iNiedzielski et al.]|2009b[ [Bowler et al.l 
I2010D . Also intriguing is the growin g population of 
brown dwarfs orbiting these stars Ce.g. iLiu et al.l 120081 : 
INiedzielski et al.ll2009bl : lQuirrenbachll2011| ). 

The most dramatic evidence for the dynamical changes 
in planetary systems around evolved stars has been pro- 
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vided by the observed absence of planets in orbits with 
semi- major axes smaller than 0.78 AU (jJohnson et alj 
|2007[) . This particular property of planets around red 
giants detected so far seems to be an imprint of stellar 
evolution on planetary systems architectures. Namely, 
the planets at small orbital separations will tidally inter- 
act with their host stars already during their MS lifetime, 
and some of them will likely end up sp iraling into their 
envelopes (see e.g. iJackson et al]|2009() . However, the 
more distant planets, which have survived the MS evo- 
lution, may face the similar fate, when their parent stars 
ascend the RGB and their tidal influe nce extends a large 
distance beyond the stellar radius (jVillaver and Liviol 
[20091) . 

A deficit of massive pla nets with short orbita l peri- 
ods was first noted by Zuc ker and Mazehl ()2002| ). The 
planet shortage in the separation range of 0.08-0.6 AU 
was explained by Udrv et al. (2003 ) as a signature of 
a transition region between two categories of planets 
that have experienced different migration scenarios. The 
gap in semi-major axis distribution fo r planets around 
M>1.2Mq stars was considered by iBurkert and Idal 
(|2007t ) who showed that the deficiency of close-in com- 
panions to giants may refiect the dearth of 0.1-0.6 AU 
companions to M^, >1.2M0 F-dwarfs caused by a shorter 
viscous diffus ion timescale of their protoplanetary disks. 
iCurrid ()2009D has shown through a series of Monte Carlo 
simulations that a stellar-mass-dependent gas disk life- 
time can explain a shortage of planets with a < 0.5 AU 
around massive stars. 

Paucity of planets orbiting clump giants within q = 
a{\ — e) = -68 A U periastron distance was first noted by 
iSato et all ()2008D . These authors singled out the planet 
engulfment at th e tip of the RGB as the m echanism re- 
sponsible for it. IVillaver and Liviol (|2009( ) have shown 
that the evolution of the star alone can quantitatively 
explain the observed absence of close-in planets around 
evolved stars, when tidal interaction is taken into ac- 
count. They have also provided the minimum orbital 
radii inside of which planets would be engulfed by stars 
of various masses at the end of the RGB evolution. These 
authors also found out that the more massive the planet, 
the earlier it w i ll be captured by the stellar envelope. 
iNordhaus et al.l (|2010[) studied in detail 1 M0 star evo- 
lution from the ZAMS to the end of the post-MS phase 
and discussed the infiuence of tidal and mass-loss laws on 
the peri od gaps in orbits of 1, 10 and 100 Mj companions. 
Finally, iKunitomo et al.l (|2011l ) found that minimum or- 
bital radii are sensitive to stellar mass in the range 1.7- 
2.1 M0 and that all the known planets indeed orbit GK 
clump giants beyond these survival limits, in ageeement 
with the planet-engulfment hypothesis. These authors 
note, however, that planets around stars more massive 
than 2.1 Mq have semi-major axes that are significantly 
larger than the minimum value, and suggest that engulf- 
ment may not be the main reason for the observed lack 
of short-period giant planets. 

In this paper, we present the discovery of a planet that 
orbits its parent star just inside the currently established 
~0.6 AU "zone of avoidance", and another one that has 
the orbital radius that is significantly smaller than 1 AU. 
We also describe in some detail the method of data analy- 
sis that we have developed over the duration of the FTPS 
project and the RV measurements of two test stars with 



known planets. 

In Section 2 we describe our observing setup and give 
parameters of the stars discussed here. In Section 3 we 
present the method of data analysis and discussion of 
the long term stability of the HET/HRS data acquisi- 
tion. The analysis of RV measurements of BD-l-15 2940 
and HD 233604 is given in Section 4, followed by the ac- 
companying analysis of rotation and stellar activity indi- 
cators in Section 5. Finally, our results are summarized 
and further discussed in Section 6. 

2. OBSERVATIONS AND STELLAR PARAMETERS 

Observations presented in this pa per were made with 
the Hobby-Eberly Telescope (HET) (JRamsev et al.|[T998l) 
equipped with the High Resolution Spectrogra ph (HRS) 
jTull 1998) in the queue scheduled mode (Shct rone et al.l 
2007). The configuration and observing procedure em- 
plo yed in our prog r am w ere identical to those described 
by iCochran et al.r(|2004[ ). The spectrograph, fed with 
the 2 arcsec fiber, was used in the R — 60,000 resolution 
mode with a gas cell (I2) inserted into the optical path. 

The spectra consisted of 46 echelle orders recorded on 
the "blue" CCD chip (407.6-592 nm) and 24 orders on 
the "red" one (602-783.8 nm). The spectral data used 
for RV measurements were extracted from the 17 or- 
ders, which cover the 505 to 592 nm range of the I2 cell 
spectrum. The spectra were reduced with the dedicated 
IRAF/Python pipeline. 

2.L BD+15 2940 

BD+15 2940 (HIC 78407) i s a ^r = 9-185 ± .024, 
B - V = 1.006 ± 0.036 (|van LeeuwenI [2001 KO 
(jCannon and Pickering! [l924[ ) star in Serpens. Its Hip- 
parcos parallax tt = 1.71 ± 1.33 points to a distance of 
585 pc. The star belongs to the FTPS Red Clump Giants 
sample (Zielihski et al. 2012). We collected 38 spectra of 
BD-K15 2940 over 2362 days between 2005 March 23 and 

2011 September 10 (MJD 53452-55814). The signal-to- 
noise ratio per resolution element (S/N) in the spectra 
was 120-330 at 594 nm in 316-919 seconds of integration, 
depending on the atmospheric conditions. 

2.2. HD 233604 

HD 233604 (BD-l-54 1280) is Vt = 10.406 ± 
0.049 (IHog et al.l [200(1. B - V ^ 1.007 ± 0.115 
(iPerrvman and ESAI Il997f ). K5 (jCannon and Pickering 
Il924( ) star in Ursa Major. It belongs to the FTPS Red 
Clump Giant sample and it was observed at 33 epochs 
over the period of 2956 days from 2004 January 11 to 

2012 February 14 (MJD 53015-55971). The S/N in the 
spectra was 110-300 at 594 nm in 762-2100 seconds of 
integration, depending on the atmospheric conditions. 

2.3. HD 209458 

HD 209458 (BD-hl8 4917, V376 Peg) a Vt = J-703 
± 0.012, B -V = 0.594 ± O 015 (Ivan LeeuwenI [2001 
GO dwarf fMontes et al ]j2001f) . is known (jHenrv et all 
2000; Mazeh et al. 2003 as a host of a transiting 
(Cha rbonneau et al.l 1200 0). to sin z = 0.689 ± 0.024 M.j 
plane t at a = 0.04723 ± 0.00079 AU orbi t ([Butler et ai.l 
[200^ (P = 3.52474859 ± 3.8x10^^ d ([Knutson et all 
120071)) with RV semi-amplitude oi K = 84.67 ± 0.7 ms'^ 
(, Torres et al..,200& ). This star is one of the FTPS test 
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stars, which are known hosts to planetary systems that 
we have been monitoring as part of our long-term ob- 
serving program. We have observed HD 209458 between 
2004 July 11 and 2011 November 25 over 2693 days or 
7.4 years (MJD 53197-55890) and collected 29 epochs of 
RV measurements. The exposure times varied between 
94 and 608 seconds depending on weather conditions, 
typically achieving the S/N of about 200. 

2.4. HD 88133 

HD 88133 (BD-H18 23 26), a Vt =.8.094 ± 0.013, 
B -V = 0.810 ± 0.015 (|van Leeuwenl [200l G5 dwarf 
(jCannon and Pickering|[l924l ) . is another FTPS test star, 
a well-studied short-period host of a msini = 0.299 ± 
0.027 Mj planet in a = 0.04717 ± 0.00079 AU (P = 
3.41587 ± 0.00059 d) with the RV semi-amplitude of 
K = 36.1 ± 3 ms-i (jFischer et all l2005l : iButler et al.l 
l2006f ). It was occasionally observed between 2006 Jan- 
uary 23 and 2012 February 20 (MJD 53758-55977). 23 
epochs of observations were gathered with S/N of 110- 
330 achieved in 141-800 seconds exposures, depending on 
actual weather conditions. 

2.5. HD 166435 

HD 166435 (BD-l-29 31 90), a Vr = 6.901 ± 0.010, 
B - V = 0.633 ±0-006 (|van LeeuwenI [200l Gl sub- 
dwarf ([White et al.ll2007f) , has been known to exhibit a 
pronounced stellar activity due to spots that manifests 
themselves as the measured RV variations induced by the 
line bisector variability of an amplitude of ^135 ms~^ 
(jCJueloz et al. 2001). As one of the FTPS test stars, 
it was observed between 2011 February 19 and 2011 
September 22 (MJD 55611-55826) at 23 epochs. The 
observations resulted in S/N of 290-540 obtained with 
70-700 seconds exposures, depending on actual weather 
conditions. 

2.6. Properties of the stars 

The atmospheric parameters {Tefj, log 5, [Fe/H], 
VmAcrn) of BD+15 29 40 and HD 233604 were taken from 
iZieliiiski et al.l (|2012| ) , who have es timate d thei r values 
using the method of iTakeda et"aLl (|200l [2001 . With 
these values and the luminosities estimated from avail- 
able data, stellar masses and ages were derived by fit- 
ting the ensemble of parameters characterizing the star 
(logi*/L0, log Tpf f, logq, and metal licity) to the evolu- 
tionary tracks of iGirardi et al.l (J2000) . The stellar radii 
were estimated using the derived masses and logq from 
the sp ectroscopic analysis as described in lZielihski et al.l 
(|2012D . The stellar rotation velocities of BD-fl5 2940 
and HD 233604 were taken from Nowak et al. (2013, 
in preparation), who have estimated their values using 
the cr oss-c orrelation technique and cali brations of Fekel! 
(|1997[) and lHekker and Mel(§nde3 (|2007[ ). The stellar pa- 
rameters are summarized in Table [TJ 

For HD 209458 and HD 8 8133 we adpoted the stel - 
lar parameters determined in iValenti and Fisched (|2005[ ) 
and lFischer et a l. (2005), respectively. 

2.6.1. Li abundance 

Lithium abundance analysis was performed 
with the Spec troscopy Made Easy package (SME, 
IValenti and Piskunov.,1996.) . which requires as an input 



a set of lines from the Vi enna Atomic Line Database 
(VALD, iKupka et al.l 119991) identified in the spectrum, 
the stellar data (including RVs), and the instrumental 
profile. Instead of the [M/H ] ratio, also ne eded by the 
SME, the [Fe/H] values from iZielihski et al.l (|2012| ) were 
used as the best approximation. The 6695-6725 A range, 
which covers both the ^Li line at 6708 A and a set of Al, 
Ti, Si and Ca spectral lines, was modeled with the SME 
using the RV, rotation velocity, and macroturbulence 
velocity as parameters. The uncertainty of the lithium 
abundance estimates was derived from the rms scatter 
of the SME fits to all the spectra available for a given 
star. In addition, the analysis of lithium abundances 
requires a correction for non-LTE effects, especially 
in the case of the lithium line at 6708 A considered 
here, which is affect ed by several non-LTE processes 
(jCarlsson et al.l [l99^ and references therein). We have 
applied the non-LTE correction from iLind et aD ()2009l ) 
to account for that effect. We have derived the values 
of A(Li) of 1.595 ± 0.093, 3.104 ± 0.056, and 1.275 
± 0.188 for HD 88133, HD 209458 and HD 166435, 
respectively. Given their evolutionary status, the Li 
content of these stars is not abnormal. For BD-l-15 2940 
we have measured A(Li) = 0.247 ± 0.132, while for 
HD 233604 A(Li) = 1.400 ± 0.042, which makes it a 
Li-reach giant. 

3. DATA ANALYSIS 
3.1. Radial velocity and line bisector measurements 

As the HRS is a general purpose spectrograph, which is 
neither temperature nor pressure controlled, calibration 
of RV measurements with this instrument is best accom- 
plished with the I2 cell technique. We have improved the 
efficiency of t he entire data analy sis process b y combinin; 
this method (|Marcv and Butleij] T992: Butler e t al.lll99' 
with that of cross-correlation (jQueloz. .1995; .Pepe et al.l 
I2OO2D for the purpose of the respective RV, line bisec- 
tor span (BS) and line bisector curvature (BC) measure- 
ments in the final version of the ALICE code. 

A detailed implementation of the above approach is 
described in iNowak (2013 ). Briefiy, we used the I2 spec- 
trum imprinted on the HRS flat-field spectrum to cor- 
rect the initial ThAr lamp-derived wavelength scale and 
to determine the instrumental profile (IP) modeled as 
a sum of five Gaussian profiles. The wavelength scale 
was assumed to be linear over the 96-pixel segments that 
the spectra were divided into. The high S/N stellar tem- 
plate spe ctrum was then s moothed using optimal Wiener 
filtering (iRabiner and Go ld 1975), deconvolved with the 
model IP using lJanssonI ( 1984il method, resampled to the 
corrected wavelength scale, and used as a model for RV 
determination from a remaining shift between the model 
and the stellar spectrum observed at a given epoch. The 
modeling process employed the non-l inear least square s 
Levenberg-Marquardt (L-M) method ([Press et al.lll992J) . 
The RV for each epoch was derived as a mean value of 391 
independent measurements from the 17 usable echelle or- 
ders, after rejecting those with outlying radial velocities 
(more than 3a off the mean value) and those obtained 
for segments without stellar spectral lines, identified in 
course of cleaning the combined stellar and iodine spectra 
from the I2 lines. The RV uncertainty for each epoch was 
calculated as the error of the mean value (= rms/yHT^, 
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where Useg is the number of actually used segments). 

A computation of the cross-correlation functions 
(CCF) for the purpose of BS/BC determination was car- 
ried out as part of the above RV measurement procedure. 
To this end, the combined stellar and iodine spectra were 
first cleaned of the I2 lines by dividing them by the cor- 
responding iodine spectra imprinted in a flat-field, and 
then cross-correlated with a binary CCF mask. Typically 
the cleaned spectra and stellar template agree within 1%, 
consistent with the S/N ratio of 100. The mask wjiscon- 
structed from a synthetic K2 star spectrum of Kiirucz: 
([1993) Atlas 9 and contained about 300 lines. 

The CCFs were computed in all the 17 HRS orders 
used for RV measurements over a wide, ± 30 kms~^ 
range, after correcting the spectra for their absolute 
radial veloc ity and shiftin g them to the Solar System 
barycenter (|Stumpfa[T980l) . Finally, the 17 CCFs were 
added together in the wavelength scale to form the inte- 
grated CCF used in further analysis. The absolute RVs 
at a given epoch of observations were measured from 
the position of a maximum of the CCF, while the CCF 
line bisectors and curvatures, defined in the usual way, 
were calculated using 3 reference ranges of (0.1 - 0.25), 
(0.375 - 0.525) and (0.65 - 0.8) of the central CCF in- 
tensity. The uncertainties of BS measurements were es- 
timate d with the formulae of iMartinez Fiorenzano et al.l 
((20051) . 

3.2. Modeling of radial velocity variations 

We fitted Kcplcrian orbits to the observed RV varia 
tions using a hybrid approach (Gozdziewski et al. 200 



Gozdziewski and Migaszewskil l2006t iGozdziewski et ah. 
2007() . which combines the global search for a plausible 



range of orbital parameters with the aid of the genetic 
algorithm (GA), followed by a non- linear least squares 
or a downhill simplex fit to data to quickly converge to 
the preliminarily constrained x^ minimum. 

In our implementation of this approach, once a periodic 
signal in a RV time series was identified in its Lomb - 
Scargle (L-S) periodogram ((Lombl 119761 : IScargld [19821 ) . 
we searched for possible orbital solutions over a wide 
range of pa rameters with the PIKAIA based GA code 
([Charbonne au 1995). The GA semi-global search usu- 
ally identified a narrow parameter range of the search 
space, which was then explored using the L-M algorithm 
to locate the best-fit Keplerian solution. I n the case of 
a single planet the RV scrambling method (I Butler et al.l 
12004) was used to assess the false alert probability (FAP) 
of the final orbital solution. If a more complicated system 
was indicated by a presence of two or more periodicities 
in the RV data, the scrambled residua were used to quan- 
tify the goodness of the final solution. Uncertainties of 
the best-fit orbital parameters were estimated from ei- 
ther the covariance matrix of the L-M fits or from the 
scrambled residuals ([Marcv et al.ll2005l) . 

3.3. Monitoring the long term stability of the 
HET/HRS data acquisition 

Along with the target stars, we have been monitoring 
three stars of known RV behavior, HD 209458, HD 88133, 
and HD 166435, to track any possible changes in the data 
acquisition system that could affect the precision of RV 
measurements . 



The RV and BS measurements for HD 209458 are listed 
in Table[2[ The long-term RV uncertainty of 12 ms~^ for 
this star was dominated by its high Te//, which resulted 
in a correspondingly low number of spectral lines that 
could be used in RV measurements. The uncertainty 
of the BS estimates amounted to 38.9 ms"""^. The L- 
S periodogram of the observed RV variations revealed 
a clear, single period of 3.52 days and no trace of any 
long-term trend. The best fit Keplerian model of these 
data is shown in Figure [T[ and the corresponding orbital 
parameters are given in Table [S[ Within uncer tainties, 
our results are identical to those of IButler et al.l ((2006). 

The RV and BS measurements for HD 88133 are listed 
in Table [3| In the case of this star, the 6.2 ms~^ RV un- 
certainty was much lower than for HD 209458, and the re- 
spective uncertainty of BS measurements was 16.5 ms~^. 
The L-S periodogram of the RV time series exhibited a 
highly significant signal at a 3.41-day period, and, as in 
the case of HD 209458, no long-period trend. The best- fit 
Keplerian orbit, folded at the orbital period is presented 
in Figure [2| and the orbital parameters are listed in Ta- 
ble [3 Once again, within uncertainties, our results are 
identical to those of Butler et al. (2006). 

Finally, despite the fact that the 51.5 ms^^ uncer- 
tainty of the BS measurements for HD 166435 (Table [H) 
was much larger than the respective RV uncertainty of 
12.7 ms~^, the correlation between the two quantities 
was obvious (r = 0.63, while the critical value of the Pear- 
son correlation coefficient at the confidence level of 0.01 
is r2i,o.oi = 0.52, Figure [3[), in ag reement with the results 
reported bv lQueloz et al.l ()2001f ) for this star. This con- 
firms the reliability of our BS measurements. Unfortu- 
nately, given the poor phase coverage of our observations 
we were not in position to confirm the periodic character 
of the RV and BS variations that mimic a planet. 

Overall, the results from RV and BS monitoring of 
the three stars discussed above demonstrate a satisfying 
long-term stability of our data acquisition scheme. 

4. RESULTS 
4.L BD+15 2940 

RV measurements of BD-f 15 2940 made at 38 epochs 
over a period of almost 6.5 years are listed in Table ^ 
and shown in Figure [1[ The estimated RV uncertainty 
for this star was ajiv = 7 ms~^ with the expected am- 
plitude of solar-type oscilla tions of Kqsc = 22 ms~^ 
([Kieldsen and Bedding||199"5l) . The star exhibits RV vari- 
ations with an amplitude of 65 ms~^ which is ~3 times 
the expected RV uncertainty. 

The L-S periodogram of the RV variations for this star 
revealed a single period of 137.5 days (Figure[n|a) and no 
signature of any low frequency component. The best-fit 
Keplerian orbit is shown in Figure [H and its parameters 
are listed in Table [5[ When an additional scatter a+ of 
18 ms~^ determined from this fit was added in quadra- 
ture to the original RV uncertainties, it resulted in a xt 
value of 0.989. Note that this increase in the long-term 
RV uncertainty is very close to the above Kosc estimate 



and the actual jitter derived as aj — \/ "o-c ~ " RV- 

The FAP<0.001 for this orbital solution was estimated 
by performing 1000 Keplerian fits to the scrambled data 
points as shown in Figure [8[ 
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Given the estimated stehar mass of 1.1 M© and the or- 
bital period of P = 137.5 days, the semi-major axis of the 
orbit of the planet is 0.54 AU, and the RV semi-amplitude 
of K — 43 ms^^ gives its minimum mass of 1.1 Mj. 
The moderately eccentric orbit of e = 0.26 ± 0.10 makes 
the planet's distance to its star vary within 0.4-0.68 AU. 
Therefore, the BD+15 2940 planet is the closest one to 
a giant star detected so far by means of the Doppler ve- 
locity method, and the second least-massive one of all 
planets orb iting giants within 1 AU after BD+48 738 b 
IjGettel et al.,,2012b ). 

4.2. HD 233604 

Radial velocities of HD 233604 (BD-l-54 1280) mea- 
sured at 33 epochs over the period of 2956 days (over 
8 years) are listed in Table [7] and presented in Figure [5l 
The estimated RV uncertainty for this star was 9.0 m s~^ , 
which is exactly the same as the estimated Kqsc value. 
The measured RV amplitude was 213 m s~^ or ~ 17 times 
the expected RV uncertainty. 

The dominant feature in the periodogram of these RV 
data is the peak at a frequency corresponding to a 192- 
day period. The corresponding best- fit Keplerian orbit, 
presented in Table[8]and in Figure[5l gives the value of xl 
of 2.76. An additional scatter of cr+=9 ms~^ was added 
to the original errors to account for an unresolved RV 
variability that has apparently degraded the precision of 
RV measurements for this star. As in the previous case, 
the PAP of <0.001 (Pigure[5]) was calculated for this best- 
fit model. The orbital solution for HD 233604, together 
with its estimated mass of Mi, — 1.5 Mq gives a msini 
= 6.6 Mj companion in a 0.75 AU, nearly circular orbit 
with e = 0.05 ± 0.03. 

5. STELLAR ACTIVITY ANALYSIS: LINE BISECTORS AND 
STELLAR PHOTOMETRY. 

In order to verify that the observed RV periodicities 
are indeed caused by the Keplerian motion, we have thor- 
oughly examined the existing photometry data in search 
for any periodic light variations, and performed a com- 
plete analysis of line bisectors and curvatures for stars 
discussed in this paper. 

5.1. BD+15 2940 

With the estimated rotation velocity of ?;rot sin z^ = 
2.3 ± 0.8 kms~i and the radius of Ri,/RQ = 14.7 ± 
2.8 the expected upper limit to the rotation period of 
BD-f 15 2940 is 317 days, but it may range between 196 
and 438 days. Therefore, it is important to consider pos- 
sible sources of stellar activity that may result in the ob- 
served periodic RV signal. Such processes include stellar 
pulsations and a spot rotating with a star. 

The p-mode oscillations can be excluded, because their 
expected amplitude and pe riod are Kpsc = 22 nis~^ an d 
Pose = 0.68 d, respectively (jKieldsen and BeddindllQQSl ). 
We do not resolve such the fast oscillation with our sam- 
pling that is focused on long-term RV variations. Con- 
sequently, this effect only contribute a RV-jitter to our 
measurements . 

The mean values of the BS and the BC for this star 
are -5 ± 23 ms~^, and -20 ± 32 ms~^, respectively, with 
the corresponding uncertainties of 18 m s~^ and 24 m s~^ . 
The BS and BC amplitudes are 45 ms~^ and 75 ms~^. 



Neither BS (Figure Ulc) nor BC show any periodic signal 
above the f„ = 0.01 level. Pearson's correlation coeffi- 
cients for RV and BS and BC are 0.12 and 0.17, respec- 
tively, with the critical value at the confidence level of 
0.01 is r36,o.oi = 0.41. 

Two archival photometry datasets are available for 
BD-l-15 2940. The exi sting Hipparcos photometry 
(jPerrvman and ESA||1997[ ) includes 65 epochs that uni- 
formly cover the 1159 days (3.2 years) or ~^8.5 RV periods 
between MJD 47904-49063, over 12 years earlier than our 
RV data. No periodic signal is present in the se data, ei- 
ther. The 362 epochs of ASAS photometry (jPoimanskil 
[T997f) collected between MJD 52685-55089 cover ~6.5 
years (2404 days) or ~12 RV periods. These measure- 
ments, clumped in 7 groups, are contemporaneous with 
^70% of our RV data. Again, no trace of a periodicity 
is present in these data (see Figure |6]d). 

Since the ASAS photometric uncertainty oi a — 0.015 
mag is typical given the observed brightness of the star, it 
is unlikely that it may be a noise caused by a stellar spot. 
However, if we made such an assumption and postulated 
a rotating spot that covers 1.5% of the stellar surface, 
the expected semi-amplitudes of RV and BS variations 
would be 61-95 ms~^ and 12-1 9 ms~^, respectively, de- 
pending on t he assumed m odel jSaar and Donahudl 19971 : 
lHatzesll200^ : iDesort et al.l l2007[). Although these values 
are in a reasonable agreement with our RV and BS data, 
there is no periodicity in the ASAS data gathered over 
^12 RV periods, and no correlation between our RV and 
BS measurements. Adding to these results a moderate 
eccentricity of the assumed Keplerian orbit, we can prac- 
tically rule out a possibility that the observed RV period- 
icity is caused by any of the stellar phenomena discussed 
above and is most likely the effect of an orbital motion 
of a substellar mass companion to the star. 

5.2. HD 233604 

The estimated upper limit to the rotation period of 
HD 233604 is 282 days and it is fairly close to the ob- 
served period of the RV variations of this star. Given 
the uncertainty in stellar radius and the precision of ro- 
tation velocity determination, the rotation period may 
range from 181 to 383 days. Obviously, it is critically 
important to examine the available data for a possible 
alternative explanation of the periodic RV signal. 

Due to the length of the period of the observed RV 
variations and their high amplitude we can confidently 
reject p-mode oscillations as a cause of the RV periodic- 
ity. The most likely alternative is a spot on the stellar 
surface moving with the star at its rotation period. Such 
a spot should be detectable in both the photometric and 
the line profile variations. 

The amplitudes and mean values of the BS and BC 
measurements for this star are 65 ms~^ and 9.6 ± 
30.1 ms^^, and 52 ms~^ and -9.7 ± 32.3 ms~^, respec- 
tively, and the corresponding uncertainties are 22 and 30 
ms~^. Neither BS (Figure [7] c) nor BC variations show 
any evidence for a periodicity in the computed LS pe- 
riodograms. Similarly, no correlations between the RV 
and the BS and BC variability were found. The respec- 
tive Pearson's correlation coefficients are r = 0.34 for the 
BS and r = 0.20 for the BC (with critical value of 7'36,o.oi 
= 0.44 for the confidence level of 0.01). 

In the case of this star, long-term photometric time se- 
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ries are availabl e from the NSVS ([Wozniak et al.l 120041) 
and the WASP (IPollacco et alJl2006n . SWASP photom- 
etry consists of 2767 measurements between MJD 54189 
and 54575, which are contemporaneous with our spectro- 
scopic observations but cover only a small part of their 
time baseline. In fact, most of the available epochs of 
photometry fall within 175 days between MJD 54400 
through MJD 54575. The NSVS photometry is of a 
better quality and consists of 309 measurements made 
between MJD 51272 and MJD 51633, which amounts 
to 361 days of data collected about 4 years earlier than 
our HET observations. These data are clustered in two 
groups, each covering a fraction of the orbital period. 

The LS analysis reveals no periodic signal in SWASP 
data, but four significant periodicities are found in the 
NSVS photometry with f„ < 0.001. Two of them that are 
close to a 1-day period are most likely artifacts of the par- 
ticular cadence of the RV observations. Two other two at 
28.1 and 27.7 days are most likely caused by the residual 
lunar illumination of the field that was not removed in 
the photometric data reduction process. 

As in the case of BD-l-15 2940, we can make a very 
conservative assumption that the SWASP photometric 
uncertainty of cr = 0.018 mag is caused by a spot, 
covering 1.8% of the stellar surface. This gives the 
expected semi-amplitude of RV and BS variations of 
70-111 ms~^ and 13-20 ms~^, respectively, depending 
on the assumed mo del (Saar and Donahue 1997; Hatzcs 
12002"; Dcsort e t aP [20071 ). The expected RV variations 
estimated under such an unfavorable assumption are 1.6- 
2.5 times smaller than those observed. 

Altogether, the absence of the periodicity seen in the 
RV data in the available photometric data, and no cor- 
relation between EC, BS and RV variations make it un- 
likely that they are induced by stellar activity. An ad- 
ditional evidence supporting the Keplerian origin of the 
observed periodic RV variations is provided by the fact 
that they are much larger than expected under the as- 
sumption that the photometric data scatter is dominated 
by a rotating spot on the stellar surface. 

6. DISCUSSION 

In this paper we have presented the detection of two 
planetary-mass companions to two red giants in the solar 
to intermediate- mass range. HD 233604 hosts a to sin i 
= 6.6 Mj planet in a 0.75 AU or ~15i?* orbit. The 
planet's closest approach to the star at the periastron 
distance of q is 0.71 AU or --14 stellar radii. BD-K15 2940 
harbors a msini = 1.1 Mj planet in a 0.54 AU or only 
~7.9 i?* orbit which approaches the star as close a.s q = 
0.4 AU or ^--5.8 i?*. Both giants belong to the Red Giant 
Clump and, given a large uncertainty in their luminosity 
determinations, their exact evolutionary stages and ages 
are not known. They are most likely evolving upwards 
the red giant branch toward the helium-flash, but their 
location on the horizontal branch cannot be excluded. 

The case of the M^ = 1.5 M© star, HD 233604, is of 
special interest, as this giant exhibits a relatively high Li 
abundance. This suggests that t he s tar is still before the 
first dredge-up process (Iben 19 67b[ ). Alternatively, but 
less likely, it may be another case of a re cent engulfment 
episo de similar to that of BD-l-48 740 (jAdamow et al.l 
l2012f) . The more likely pre-dredge-up scenario would 
mean that the HD 233604 planetary system was not yet 



affected by a substantial stellar radius growth before the 
helium flash and the subsequent severe mass-loss. The 
low rotation velocity of HD 233604 additionally supports 
this scenario. A more detailed, ongoing chemical com- 
position study will allow constraining the evolutionary 
status of this star much better. 

The status of the M^, = 1.1 Mq star, BD-H15 2940, is 
less certain. Its poorly determined but definitely rela- 
tively high luminosity suggests that the star is a rather 
evolved. He-burning object already past the episode of 
stellar radius expansion. 

A useful way to illustrate the occurrence of tight-orbit 
planets around various types of stars is to plot their semi- 
major-axes against the values of the stel lar logg. Such a 
plot o f the current catalog of 778 planets (" Schneider et al.l 
120111 ) is shown in Figure [HI Clearly, the most recent data 
show that there are no planets within the borderline of 
log{a[AU]) = -0.633 logg-f 0.975 (set by HD 96127 b and 
Kepler- 10 b positions in that graph). Th e only exception 
is th e intriguing system of HIP 13044 (jSetiawan et al.l 

[2oToh . 

Planets over the whole range of orbital separations, 
all the way down to Hot Jupiters, can be found around 
MS stars with 4< logg <5. Also, sub-giant (logg ^ 
3.5) companions still cover a wide range of orbital radii, 
with the curre nt lower limit of 0.08 10 AU in the case 
of HD 102956 (j Johnson et al.ll2010bf) . but, interestingly, 
there are very few planets inside 1 AU of these stars. In 
particular none is known within 0.08 AU< a <0.76 AU 
in ac cordance with the predictions of (jBurkert and Idal 
12001 . 

Planets in orbits around more evolved stars with even 
lower logg ^2-3 are indeed found outside the pre- 
dicted borderline, as illustrated by the 0.54 AU orbit 
of the BD-l-15 2940 planet discussed here. However, 
one notes that planets around giants are generally found 
in tighter orbits than those around sub-giants (except 
HD 102956 b). This suggests that, in addition to the 
absence of 0.1-0.6 AU companions to Mi, >1.2 Mq F 
dwarfs, another effect must influence the orbits of plan- 
ets around those evolved stars. This effect has been 
identified as the tidal interaction between the star and 
the planet, whic h may cause an accelerated orbital de- 
cay (JRasio et al.lll996HVillaver and Livi6l l2009D . Accord- 
ing to lVillaver and Livid (|2009D . a 2 M© star, a typical 
object in the observed population within the logg ^2- 
3 range, should have already ingested 3-5 Mj or more 
massive planets in orbits up to 0.2 AU before reaching 
\ogLi,/LQ = 2 (or logg ~2.3) on the Red Giant Branch. 

Both HD 233604 b and BD-hl5 2940 b fall into that 
part of Figure [9] as well. However, HD 233604 b in 
0.75 AU or of g = 0.71 AU orbit around its M* = 1.5 Mg 
host is located well within th e minim um orbital radius as 
define d in Villaver and Livid ()2009f ) and lKunitomo et al.l 
(|2011l ) for a Jupiter mass planet. With the estimated 
M^ = 1.1 M0 BD-hl5 2940 planet in 0.54 AU or only q 
= 0.4 AU i s located even closer, wi thin the star's RGB 
tip radius (jVillaver and LivTol 1200 9^ . Existence of such 
close-in planets suggests that the influence of star-planet 
tidal interactions in late stages of stellar evolution may 
be overestimated. The relatively large eccentricity, e = 
0.26 ± 0.1, of the BD-hl5 2940 planet also deserves more 
attention as one might expect its orbit to be already cir- 
cularized at the distance of only 5.8 stellar radii from 
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the star. On the other hand, it may represent a reminis- 
cence of a recent plan et-planet interaction in the system 
(jVovatzis et al.ll2013[ ). 

Even more evolved hosts of planetary systems, with 
log 5 < 2, form a very heterogenous group with masses 
ranging from around 1 Mq (42 Dra) up to Mi, = 
4.5 Mq (HD 13189). As the evolution of planetary 
systems affected by tidal interactions and mass loss is 
very much stellar mass-dependent, it is difficult to dis- 
cuss them together. These stars also exhibit substantial 
amounts of variability characterized by post-fit RV rms 
values easily reaching 55-60 ms~^ (jHatzes et al.l 1200^ 
iNiedzielski et al.ll2007 ) which makes studying them even 



more difficult. 

Although the population of planetary systems detected 
around evolved stars has been steadily growing, it is still 
too meager to achieve satisfying levels of understand- 
ing of the star-planetary system interaction driven by 
the post-MS stellar evolution. An additional, impor- 
tant complication stems from the fact that luminosities 
of these stars are typically poorly determined, which ob- 
viously makes constraining their evolutionary status very 
difficult. In any case, continuing studies of such systems 
provide important contributions to the process of creat- 
ing a general picture of the evolution of planetary systems 
under the influence of their evolving parent stars. 
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TABLE 1 
Stellar parameters of BD+15 2940 and HD 233604. 



Parameter 



BD+15 2940 



HD 233604 



Vt (mag) 

{B - V) (mag) 

Spectral type 

TV (mas) 

T.ff (K) 

logs 

(Fe/H) 

logL*/L0 

M*/Mq 

R*/Rq 

Age (Gyr) 

Vrotsmi-t (kms^^) 

Prot (sin u)^-*- (days) 

Koaa (ms~l) 

Pose (days) 



9.185 ± 0.024 
1.006 ± 0.036 
KO 

1.71 ± 1.33 
4796 ± 117 
2.8 ± 0.45 
0.28 ± 0.07 
2.01 ± 0.75 
1.1 ± 0.2 
14.7 ± 2.8 
7.1 ± 4.2 
2.3 ± 0.8 
317 ± 121 

99+137 

0.68+-^l 



10.406 ± 0.049 
1.007 ± 0.115 
K5 

4791 ± 45 
2.55 ± 0.18 
-0.36 ± 0.04 
1.75 ± 0.18 
1.50 ± 0.30 
10.9 ± 0.6 
2.4 ± 1.8 
2.0 ± 0.8 
282 ± 111 
q+9 



0.27 



+0.11 
-0.07 



TABLE 2 
Relative radial velocities and bisector span of HD 209458. 



Epoch (MJD) 



RV (ms-i) 



o-Rv (ms ^) 



BS (ms--^) 



o-BS (ms ^) 



53197. 
53538. 
53544. 
53547. 
53984. 
54001. 
54014. 
54056. 
54448. 
54738. 
54756. 
54759. 
54764. 
54777. 
54780. 
54786. 
54792. 
54813. 
55031. 
55068. 
55110. 
55502. 
55504. 
55506. 
55510. 
55510. 
55513. 
55515. 
55890. 



29482 
35438 
34874 
33435 
34746 
09852 
26729 
13726 
07560 
09154 
22903 
20674 
21486 
16241 
15300 
13373 
13299 
07756 
26789 
16183 
25281 
17681 
18200 
18306 
15573 
15868 
16310 
15433 
13839 



-80.9 
-73.1 

49.9 

57.4 

65.7 
-53.7 
-34.1 

50.6 
-59.1 

18.8 

53.4 
-32.3 
100.9 

27.9 
-21.4 
-26.2 
122.2 
115.0 

63.4 

-28.8 

3.1 

-67.2 

96.2 
-60.1 
-33.5 
-35.6 
-77.0 

34.4 
-82.0 



13.1 
14.7 
13.9 
14.5 
13.3 
11.3 

9.7 
12.1 

7.3 
15.3 
11.2 
13.4 
10.5 

9.5 
12.1 
12.4 
10.9 
10.9 
15.3 
15.0 
11.8 
12.2 

9.6 

9.5 
11.2 

9.2 
12.9 
13.6 

7.7 



39.9 

123.9 

-11.5 

-19.8 

-13.6 

19.7 

3.5 

60.2 

33.6 

-0.0 

70.3 

9.9 

-18.3 

38.0 

-14.0 

82.3 

39.0 

95.2 

18.2 

-8.5 

-61.2 

-19.0 

-23.4 

9.6 

40.4 

39.7 

47.1 

-47.7 

56.8 



37.3 
46.9 
37.9 
39.0 
40.3 
33.4 
36.2 
42.9 
26.2 
52.1 
35.9 
43.9 
44.2 
39.3 
48.3 
54.2 
44.2 
36.8 
48.1 
38.1 
38.6 
41.8 
26.1 
28.4 
39.4 
36.5 
36.8 
29.6 
25.2 
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TABLE 3 
Relative radial velocities and bisector span of HD 88133. 



Epoch (MJD) 



RV (ms-i) 



o-RV (ms 1) 



BS (ms-i) 



o-Bs (ms 1) 



53758. 
53774. 
53835. 
53868. 
54090. 
54090. 
54165. 
54461. 
54465. 
54474. 
54767. 
54777. 
54782. 
54785. 
54788. 
54858. 
54864. 
55137. 
55940. 
55948. 
55956. 
55973. 
55977. 



26950 
21259 
25918 
17657 
36895 
37559 
36426 
34054 
32153 
31236 
50684 
48514 
47424 
45165 
45604 
46134 
44610 
50172 
50047 
27405 
24389 
40238 
19226 



-21.4 

33.2 

27.5 

-19.8 

-31.8 

-30.5 

-21.8 

8.9 

-32.2 

43.2 

29.3 

12.0 

-7.1 

18.8 

29.1 

-28.6 

-9.5 

-15.8 

-34.2 

-13.6 

37.6 

28.1 

12.5 



5.5 
6.0 
6.6 
6.5 
5.3 
6.6 
6.3 
5.6 
5.6 
6.2 
7.5 
6.2 
6.1 
6.5 
5.4 
8.1 
7.4 
5.1 
6.1 
4.7 
5.4 
6.4 
6.3 



-28.8 
-30.2 
-27.5 
-12.8 
-18.4 
-43.0 

15.0 
-32.0 
-13.6 
-24.9 

22.7 
-41.6 

28.8 

-5.3 

9.4 

-11.0 

-12.1 

-8.7 

2.9 

-60.2 

11.4 

-9.5 
-11.8 



14.2 
14.6 
13.4 
13.8 
17.6 
19.8 
16.4 
15.1 
13.5 
16.0 
19.4 
20.2 
19.6 
18.8 
16.7 
21.5 
19.1 
14.2 
17.0 
11.1 
10.9 
18.9 
16.9 



TABLE 4 
Relative radial velocities and bisector span of HD 166435. 



Epoch (MJD) 



RV (ms-i) 



o-Rv (ms '-) 



BS (ms-i) 



o-BS (ms '-) 



55611. 
55612. 
55613. 
55615. 
55616. 
55626. 
55627. 
55628. 
55634. 
55639. 
55640. 
55643. 
55644. 
55645. 
55801. 
55805. 
55814. 
55815. 
55816. 
55817. 
55824. 
55825. 
55826. 



51710 
50123 
51535 
51164 
50156 
49751 
47291 
49111 
45595 
44095 
44137 
41648 
42200 
40177 
24069 
21921 
18152 
18623 
18358 
18530 
15294 
14860 
13646 



-33.2 

-33.4 

10.5 

-24.1 

29.1 

-15.8 

15.2 

12.1 

51.8 

-7.3 

-15.3 

-41.3 

46.8 

30.3 

26.1 

43.0 

-28.4 

-59.8 

34.5 

-1.6 

4.7 

-24.6 

0.2 



10.8 
8.4 
15.7 
16.8 
10.2 
18.0 
10.3 
18.9 
10.3 
12.8 
10.2 
16.4 
15.2 
18.1 
11.0 
11.0 
9.0 
12.8 
11.0 
11.5 
11.3 
11.1 
11.5 



27.7 

-3.9 

23.3 

-54.3 

-47.1 

20.3 

-54.4 

-99.7 

-104.3 

10.9 

-77.0 

52.0 

-91.9 

-92.1 

-154.1 

-145.8 

-14.6 

-32.5 

-45.3 

-70.8 

-18.2 

-28.3 

-129.4 



56.5 
52.1 
52.9 
51.8 
58.9 
56.1 
53.9 
66.7 
55.1 
68.0 
60.8 
61.4 
44.8 
63.6 
48.4 
44.6 
40.2 
40.5 
43.3 
42.5 
41.2 
42.2 
39.5 
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TABLE 5 
Orbital parameters of HD 209458 b and HD 88133 



Parameter 



HD 209458 b 



HD 88133 b 



P (days) 
To (MJD) 
K (ms-l) 
e 

a; (deg) 
msini (Mj) 
a (AU) 
RVo (ms-i) 
(T-i- (ms^-'-) 

f^O- C (m_s"^) 
o-RV (ms^i) 
cFj (ms^-^) 

Nobs 



3.525 ± 0.009 

53199.3 ± 0.9 

87.7 ± 5.4 

0.007 ± 0.060 

60 ± 108 

0.714 

0.0474 

11.1 ± 3.4 

0.0 

3.13 

20.54 

11.87 

16.76 

29 



3.415 ± 0.001 

53760.7 ± 0.9 

34.2 ± 3.6 

0.09 ± 0.07 

15 ± 120 

0.284 

0.047 

-0.7 ± 

0.0 

1.127 

5.46 

6.15 

sS 5.46 

23 



1.4 



TABLE 6 
Relative radial velocities and bisector span of BD-I-15 2940. 



Epoch (MJD) 



RV (ms-i) 



o-Rv (ms 1) 



BS (ms-i) 



o-BS (ms 1) 



53452. 
53499. 
53844. 
53866. 
53875. 
53892. 
53902. 
53910. 
53955. 
53971. 
54215. 
54225. 
54248. 
54276. 
54290. 
54320. 
54515. 
54564. 
54595. 
54624. 
54692. 
55018. 
55049. 
55254. 
55280. 
55582. 
55621. 
55621. 
55674. 
55719. 
55742. 
55750. 
55763. 
55776. 
55795. 
55801. 
55808. 
55814. 



36375 
41176 
29702 
40534 
39832 
16699 
31759 
29146 
16666 
13068 
25852 
42513 
35726 
28661 
25121 
17812 
45514 
30494 
22798 
34067 
15651 
25340 
18227 
40663 
37572 
52001 
40203 
41021 
46852 
34138 
28347 
25899 
22678 
18611 
13382 
11490 
08974 
09263 



-46.8 

8.4 

-40.7 

10.7 

19.5 

37.0 

62.6 

23.5 

-20.7 

-40.3 

-32.1 

-27.4 

-29.3 

-6.9 

38.6 

37.9 

-30.1 

26.3 

47.0 

-38.2 

17.6 

22.0 

-54.8 

28.6 

27.5 

13.9 

-1.5 

-12.4 

6.9 

-4.3 

-16.6 

-37.5 

-53.9 

-68.1 

-3.2 

35.3 

60.0 

43.2 



7.8 
7.5 
5.5 
6.6 
5.6 
6.3 
5.5 
5.9 
5.2 
5.9 
6.6 
7.9 
6.1 
7.6 
7.3 
7.6 
7.6 
7.8 
7.7 
7.2 
8.2 
8.7 
6.4 
7.8 
8.0 
7.6 
7.5 
6.0 
7.5 
6.6 
6.9 
7.7 
8.5 
6.6 
6.4 
6.1 
6.7 
6.3 



15.9 

5.9 

-4.5 

-7.3 

-4.7 

25.4 

1.7 

0.3 

-45.9 

-8.1 

28.2 

-15.1 

-36.4 

42.4 

-45.6 

4.6 

-30.5 

-32.5 

9.9 

-24.2 

25.1 

-17.2 

9.2 

7.2 

34.5 

-11.6 

30.7 

16.7 

-5.5 

-29.0 

-25.8 

4.4 

-37.6 

-11.2 

-11.2 

3.0 

-17.4 

-18.1 



17.0 
17.8 
14.5 
18.5 
13.0 
16.6 
14.9 
14.4 
16.1 
14.8 
17.0 
26.1 
17.1 
19.4 
19.5 
20.9 
19.2 
23.3 
17.4 
16.3 
19.1 
17.1 
15.0 
20.8 
22.8 
20.7 
24.0 
15.0 
18.4 
15.6 
15.5 
14.6 
17.2 
16.1 
15.1 
15.4 
17.1 
17.5 
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TABLE 7 
Relative radial velocities and bisector span of HD 233604. 



Epoch (MJD) 



RV (ms-i) 



o-Rv (ms 1) 



BS (ms-i) 



o-BS (ms 1) 



53015. 
53360. 
54064. 
54136. 
54208. 
54399. 
54420. 
54500. 
54564. 
54604. 
54763. 
54782. 
54839. 
54865. 
54891. 
55173. 
55173. 
55176. 
55214. 
55223. 
55254. 
55307. 
55486. 
55487. 
55502. 
55519. 
55520. 
55645. 
55659. 
55671. 
55956. 
55964. 
55971. 



25986 
29494 
39510 
38069 
19384 
46353 
41574 
17910 
21884 
14329 
46178 
43096 
45575 
19007 
33472 
35355 
37553 
52706 
26133 
20532 
36068 
19447 
49435 
49885 
47794 
38507 
39557 
26266 
21377 
19521 
18941 
38336 
16102 



-28.9 

-231.5 

24.8 

-172.2 

142.1 

135.1 

138.1 

-193.3 

80.0 

119.2 

83.8 

193.9 

-7.2 

-134.9 

-183.9 

165.5 

159.9 

181.0 

5.7 

-28.9 

-151.2 

-102.0 

-145.8 

-155.4 

-56.7 

18.2 

35.2 

-155.3 

-190.8 

-185.3 

167.2 

152.6 

59.5 



10.6 

10.7 

7.9 

10.6 

9.2 

7.5 

8.3 

10.0 

9.0 

8.5 

9.5 

9.0 

10.1 

10.6 

11.9 

9.9 

8.6 

9.5 

9.9 

8.5 

10.3 

8.1 

7.5 

6.9 

9.7 

10.6 

8.9 

9.1 

8.2 

7.2 

6.9 

7.2 

7.5 



4.8 

-57.9 

43.7 

-2.6 

12.1 

8.3 

5.2 

-50.3 

29.2 

12.1 

-20.0 

14.1 

0.1 

15.4 

-7.2 

22.4 

1.3 

64.1 

72.5 

36.4 

-43.7 

-10.7 

13.2 

-8.3 

44.5 

35.5 

-10.1 

49.6 

26.4 

5.7 

34.3 

3.6 

-28.4 



26.2 
22.0 
19.3 
20.6 
18.5 
20.6 
24.2 
25.1 
21.4 
20.5 
26.1 
28.3 
18.1 
26.7 
25.9 
26.1 
24.0 
28.6 
25.0 
24.8 
24.1 
16.1 
21.5 
18.0 
26.9 
30.3 
17.4 
23.1 
15.6 
17.3 
17.1 
16.4 
15.9 



TABLE 8 
Orbital parameters of BD+15 2940 b and HD 233604 b. 



Parameter 



BD+15 2940 b 



HD 233604 b 



P (days) 
To (MJD) 
K (ms-l) 
e 

ui (deg) 
msinj (Mj) 
a (AU) 
RVo (ms-i) 
(T-i- (ms^-') 

(To-C (ms^i) 
CTRv (ms~^) 
cFj (ins^-^) 

Nobs 



137.48 ± 0.34 


53464 ± 16 


42.7 ± 4.4 


0.26 ± 0.10 


302 ± 33 


1.11 


0.539 


-2.5 ± 3.3 


18.0 


0.989 


17.89 


6.97 


16.48 


38 



192.00 ± 0.22 

53022 ± 77 

177.8 ± 4.3 

0.05 ± 0.03 

281 ± 43 

6.575 

0.747 

-16.9 ± 3.7 

9.0 

2.76 

19.13 

9.04 

16.86 

33 
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Fig. 1. — Radial velocities of HD 209458 plotted as a function of orbital phase for the solution detailed in Table [5] The measurements 
taken at or near the transits have been removed. 



Planetary-mass companions to the K-giants BD+15 2940 and HD 233604. 

HD 88133 

T 1 1 1 1 1 1 1 r 



15 




I 20 F 



CO 



10 ^ . 




Y -10 ^ 

o 



1 1 1 1 1 1 1 1 r 



i; S ' ii H % i 



1 ^ 



\ 



hi 



ft 



J I I I I I I I L 







0.5 




Fig. 2. — Radial velocities of HD 88133 plotted as a function of orbital phase for the solution detailed in Table [5] 
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Fig. 4. — Left: radial velocities (filled circles), the best fit Kcplcrian model detailed in Table |8] (full line) and the post-fit residuals for 
BD-f 15 2940. Right: the same as a function of orbital phase. 
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Fig. 5. — Left: radial velocities (filled circles), the best fit Kcplcrian model detailed in Table |8] (full lino) and the post-fit residuals for 
HD 233604. Right: the same as a function of orbital phase. 
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Fig. 6. — The L-S periodograms of (a) radial velocities (b) RV residua, (c) bisector span and (d) ASAS IIPoimanskilll997l) photometry of 
BD+15 2940. The levels of FAP=1.0% and 0.1% are shown. 
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Fig. 7. — The L-S periodograms of (a) radial velocities and their residua (b), (c) bisector span and (d) NSVS l lWozniak et al.l 120041 ) 
photometry of HD 233604. The levels of FAP=1.0% and 0.1% are shown. 
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Fig. 8. — Histograms of the values of (x^)^ obtained from the fits of the Keplerian models to scrambled RVs used to estimate the 
FAPs for the two stars discussed in the text. In each case, 1000 sets of scrambled RVs have been generated. Vertical arrows point to the 
respective (xi)^'^ values derived from fits to unscrambled velocities. In both cases the FAP values are less than 0.1%. 
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Fig. 9. — Semi-major axes of 778 planets as a function of log ly of their host stars. The sample of 84 planets orbiting giants and subgiants 
with R* > 3 R© (0) includes 12 planets discovered witin FTPS survey (*). BD+15 2940 b and HD 233604 b are indicated by * with a 
rim. 



